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1. INTRODUCTION 

The problem every day and around the world is to meet the energy demand that is growing more and 
more. The accelerated increase in industrial activity in countries and business investment that provides a 
lower cost of production has increased rapidly the global demand for energy. Faced with this demand, the 
industrialized countries have massively used the production of decentralized energies whether renewable or 
non-renewable sources [1, 2]. Wind energy was the first of the energy sources exploited by man. This 
inexhaustible source of energy has developed significantly in recent decades [2, 3]. 

However, the large proportion of wind turbines uses DFIG (asynchronous double-feed generators). 
This generator allows variable speed electricity generation [4]. It gives the opportunity, to better control the 
wind resources for different wind conditions [5]. In the field of variable speed, the asynchronous generator 
and thanks to the development of electronic inverters, currently ensures, in addition, a very important market 
share. In the last years, an appreciable increase in the order of generators with double feeding have been 
given, where we can give the vector control, the direct control of the couple and the adaptive control. 

Similarly, multi-scientific works on the DFIG have contributed to the controls of wind systems 
(backstepping control, Sliding mode control, MPPT control, PCHD Model control, D-FOC control etc. [1, 2, 6-10]. 
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The major disadvantage of these commands is assembled on a rigorous mathematical aspect for synthesized. In 
contrast, there are other commands that have a physical meaning such as passive control [11-14]. 

Passivity-based control is a well-established technique that has been very powerful in designing robust 
controls for physical systems. Inherited from well-known physical phenomena, notions of passivity are adapted to 
several scientific domains and are effective for the regulation of electrical, mechanical and electromechanical 
systems present in several fields of engineering, such as robotics, power electronics etc. [6, 15, 16]. 

Passivity-based control provides robust controllers that have a clear physical interpretation in terms 
of the system's interconnections with its environment. In particular, the total energy of the closed-loop system 
is the difference between the energy of the system and the energy supplied by the controller [7]. Moreover, 
since the Euler-Lagrange structure is preserved in a closed loop, the passivity control has a robust stability 
with respect to the non-modelled dissipative effects and exhibits robust performances due to its inverse 
optimality [8, 9]. Modeling in PCH (Port Controlled Hamiltonian) form and its use for non-linear system 
control has already been successfully applied in many areas of physics to describe these systems in terms of 
exchange of energy [17]. This by generalizing the notion of Dirac structure into an interconnection structure 
that makes it possible to clearly write the interaction of the system with its environment [18, 19]. 

The structure that includes power converters on the rotor part provides a high degree of 
controllability, and allows for maximum efficiency at all wind speeds. To ensure the control of MPPT 
(Maximum Power Point Tracking) we used the diagram presented in Figure 1 [20, 21]. 
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Figure 1. Block diagram with speed control 


The control scheme of the DFIG requires that the rotor speed in optimal operation, therefore, the use 
of speed sensor to ensure this magnitude implies the increase in installation complexity and cost. For this 
purpose, the Model Reference Adaptive System (MRAS) observer was used to estimate the rotor position 
from two rotor flux models (reference model, adjustable model) to obtain the rotor speed [21]. Figure 2 
presents the schematic diagram of an MRAS observer. 
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Figure 2. Principle diagram of an MRAS observer 


Nevertheless, the objective of this work shows the resolution of the three problems generally encountered 
in the control of the DFIG, then the first problem is the complexity of the mathematical equations which classically 
models the system to study, the second problem is the minimization of the costs of construction by the elimination 
of mechanical sensors, and finally we model and control the production chain in closed loop. Therefore, the 
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proposed solutions present equations that have a physical aspect (energetic) and a non linear control IDA-PBC 
(Interconnection and Damping Assignment - Passivity Based Control) without a mechanical speed sensor 
associated with a wind turbine integrated into the electrical network (Figure 3 shows the overall system studied). 
The structure of this contribution is based on the modeling and control of the system, results and discussion finally 
a conclusion and perspective. Variants and scenarios are performed under the Matlab / Simulink environment to 
evaluate the reliability and robustness of sensorless passive control. 





Figure 3. Structure of sensorless control with IDA-PBC method for the DFIG 


2. MODELING OF THE WIND SYSTEM 
2.1. Modeling of the wind turbine 
The wind power equation is defined from the following way [3]: 


1 
P, =>7p.S. v" (01) 


p : is the density of the air. 

S : is the circular surface swept by the turbine. 

v : is the wind speed (m / s). 

The aerodynamic power of the turbine is then written: 


1 
Feros = Cy. Py = Cy (A, B) = p.S.v° (02) 


With A which represents the specific speed 
AÀ = RQwind 


Vv 
The power coefficient Cp represent the aerodynamic yield of the wind turbine, it is determined as following: 


1 ” C6 
Co(A, B) = C, (C5 — CaB* — Ca — C5) e4 (03) 
Knowing that: 
1 (1 0.035 
A (z a) (04) 


The torque of the wind turbine is given by: 


A wind turbine sensorless automatic control systems, analysis, modelling and ... (Hamiani Hichem) 


48 0 ISSN: 2088-8694 


Sv’ À, 
C= Pmec _ 19-5.V"Cp(4,B) (05) 


Ww 2 w 





The fundamental equation of dynamics allows determine the evolution of the mechanical speed from the total 
mechanical torque Tmec applied to the rotor: 


dQmec 
J Ea = Tmec (06) 
J : is the total inertia that appears on the rotor of the generator. 
Jturbine 
J= = + JpFi¢ (07) 


This mechanical torque Tmec takes into account, the electromagnetic torque Tem produced by the 
generator, the Torque of viscous friction Tvis, and the torque from the gearbox Tọ. 


Tmec = Ty = Tom 7 Tyis (08) 


2.2. Modeling of the DFIG 
The model of the double-feed generator in the d-q reference can be written as below [8]: 


= d®ds 
Vas = Rslas + dt < ws Pas 


d® 
= qs 
Vas = Rel gs + dt _ ws Pas 





Var = Rrlay + EE — 0,Oqr a 
Var = Rylgr + = 0, Day 
With : 
w, = w, — PA (10) 


Das = Lslas + Mlar 
Oj = Lslgs + Mio is 
Dar = Ly lar T Mlas 
Por = Lylgr + Mlos 


The expression of the electromagnetic torque of the DFIG depending on flow and stator currents can be 
written as follows: 


M 
Tem = -P 7 (®qslar — Daslar) (12) 


With p : number of pole pairs of the DFIG. 
The active and reactive powers stator and rotor of the DFIG are written as follows [15]: 


Ps — Vaslas + Vast as 
Qs = Vastas = Vaslas (13) 
P, E Var lar + Vorlar 
Q; = Vorlar = Varlar 


2.3. IDA-PBC controller applied to DFIG generator 

From (09), it is clear that the dynamic model of the DFIG is non-linear due to the coupling between 
the speed and the electric currents. According to the vector control principle, the direct axis current Id is 
always forced to equal zero in order to orient all the linkage flux in the d axis and achieve maximum torque 
per ampere. 
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dds 
= = Vas—Rslas + Wshslgs + wsMlqr 
d®gs 
= Vas —Rslqs a WsLslas zs WsM lar 
d® gy 
= Var Relay + wrLylgr + OMI 45 
dDgy 
a = Vor —Rylor — Ly lay + wrM las 


So the mechanical equation of the rotating part of the generator is given by: 


Ja = MEhl, -Tp — Ty 

Where ; J = k z] gE | aS ia 
The state variables are : 

va [Ol OF Jol =E xp? 
With : 
xf = [T ©] : are the electric state variables. 


Xm = Jw: mechanical variable. 
The energy function written by: 





bare 1 
HC) = x L AF x 
2 2]DFIG 


With: 


mine Mh: mao 
Re 1p = |, i 


ML L-i 


The partial derivatives of energy with respect to the state variables are: 


ðH ðH 

SS Z s] =| g] 
aXe e Xe [ S a] 
ðH = J-1% OH _ 

OXm me Oxm 7 


Finally, the matrices of interconnection, damping and the matrix of the command are : 


—WsLsJg —WsMJ, 02x34 Rsh  O2x%2 02x1 
J(x) =|-@shsJg —wrLsJ2 Mhzls|, R(x) = |02x2 Rrh 02x41 
01x2 MIŤJ, 0 Oix2 Oix2 Ty 


I, 02x2 02x1 
g(x) = 02x2 l 02x1 u = [Vf Ve T,|" 
Oix2 01x2 1 
With: 


0 0 0 Va Va 
022 = lo al> O2x1 = [g] > 01x2 = [0 0], Vf = {|W = | 


J@)=J@)*, R(x) = R(x)" 2 0 


With these matrices the model PCH is written by: 
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—w@sLsJz —WsMJz 02x1 Rsl2 02x2 02x1 h 02x2 02x1] |V 
X = —w,;MJ, —W,LsJo MJ; = 02x2 Rh 02x1 VH + 03x2 I, 02x1 yT 
01x2 MI; J, 0 Qix2 O1x2 Ty Oix2 Oix2 1 I7 
(24) 
h O2x2 02x1 
Y= |02x2 h 02x1)VH (25) 
01x2 01x2 1 
2.4. Calculation of control voltages Vdr, Vqr : 
To calculate the control voltages one has to determine Ja (x) and Ra (x) 
For that we must determine J, (x), Ra(x) the controller. The closed-loop system is written by: 
fx) t+ gU = (Ja) — Ra (x))ð Ha (x) 26) 


Where : 


Ja) =J (x) + Ja (x) 
Ra(x) = R(x) + Ra (x) 
Ha(x) = H(x) + Ha (x) 
With : 
Ha(x) : The energy function of the closed-loop system, Ha (x) : The controller energy function. 
Ja(x) : The interconnection matrix of the closed-loop system, Ja (x) : The controller interconnection matrix. 


Ra(x) : The damping matrix of the closed-loop system, Ra (x) : The controller damping matrix 
It is possible to write (26) as follows: 


U(x) + Jax) — R(x) + Ra(x))OHa(x) = (Ja) — Rae) 0H (x) + g(x)u (27) 


The previous equation is equivalent to: 


(Ja) — Ra(x))0Hg(x) = —(Ja (x) — Ra(x))OH (x) + g(x)u (28) 


So the total energy desired: 





1 * = * 1 * 
Ha(x) = 5 (Xe ~ X) L “(i _ Xe) + 2J (Xm ~ Xe) (29) 
DFIG 
So: 
Tr—1 1 * T: -T —1..% 1 «2 
H (x) = Ha(x) — A(x) = -x L * Xe — —— XX tixe L *x} + x. (30) 
JDFIG 2 2] DFIG 
With: 


0H,(x) = E , Where I = [las Iqs lar Iqr |" 
Using this relationship, (28) becomes : 


(Jax) = Ra) [TE] = -Ga() - Ra) [1] + gu (31) 


The command Vr is found in line 3 and 4 of the matrix in (31). 
So: 
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02x2 02x2 02x1 02x2 02x2 02x1 
Ja(x) = [02x2 02x2 —Srm(X)] , Ra(x) = oo rh oe (32) 
01x2 Jim) 0 Oix2 Oix2 0 
Where : 
2x1 
Jm) E F to be determined. 
r : It is an additional resistance for currents to dampen transient oscillations. 
We replace the matrices J, (x),R, (x) and using (30), we find: 
T (Fara *k\T * 
Jrm (x) =M PERE CU; — Is) Joly (33) 
So: 
V = W — (w = w)(L Jai} + Jem) - Mo* Jos — 13) -rh — Tr) (34) 


Unfortunately, the order is singular in point of equilibrium. It's done possible to get rid of this 
singularity by adding a depreciation variable. Now we keep Ja(x) and H,(x) as before, and we change the 
matrix to have the form: 


02x2 02x2 O2x1 
Ra(x) = [02x2 rh O2x4 (35) 
0ix2 01x2 g(x) 
With : 
Tem—Tem e 
f(x) = “em Teme) (36) 
and : 
J ine = Cw” (37) 


When replaced in the closed-loop Hamiltonian equation, §(x) is multiplied by (w — w*), therefore 
elimination of singularity. Since the mechanical part of (27) is the only one that has been modified, 
then the expression 
of Vr in terms of J, (x) rest the same. 

Using the equilibrium equations, we obtain: 


Jim) = MJzl; (38) 


The closed-loop dynamic system is always of the form (24) with: 


—WslsJg ~—w@sMJ2 02x1 Rsl2 02x2 02x1 
Jax) = |7sM]2> —wrLsJz Mhls|; R(x) = |02x2 (Rr +N 02x1 (39) 
01x2 MI; J. 0 01x2 01x2 Tr + &(x) 


Finally, the rotor voltages of the order are written by: 

Ve = Ve — (w — @")(L,Jalr + MJols) — Mw" Jo — 15) — rh (l, — Ir) (40) 
With: 

Ve = (ws — @*) (Ly Joly + MJzls) + Rylaly (41) 


2.5. Application of the MRAS 
The equations of the chosen reference model and adaptive model are given by: 
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Reference model: 














aD gr Ly dl ds 

- = M ( Vas—Rslas _ OL, A ) (42) 
aDgr Lr dlqs 

oe ( Vos Relas -ols e ) (+3) 

Adaptive model: 

d® q+ 1 = At M 
a = Vas — T, Par a Par + 7, las (44) 
d®gr 1 = A> M 

7 = Vas = T Dor = OO jr + 7, tas (45) 


The error between the two models, function to be minimized, is given by: 
e= 0,- Ô, (46) 


The law of adaptation presents an integration in open loop (problem of offset). For the improvement of the 
estimation response a low pass filter has been proposed by many authors [10]: 


w = Kp (EqgD 9, — EqPar) + K;(€qP gr _ Egar Jdt (47) 
With: Kp, K; positive constants. 
3. RESULTS AND DISCUSSION 


Using Matlab/Simulink, the results obtained for the sensorless IDA-PBC command are presented in 
Figure 4, Figure 5, Figure 6, and Figure 7 shows the reference of the dynamic wind speed. 
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Figure 4. Results of the active power with its error, without sensor. 
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Figure 5. Results of the reactive power with its error, without sensor. 
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Figure 6. Results of the mechanical speed with its error, without sensor. 
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Figure 7. The reference of the dynamic wind speed. 


4. DISCUSSION 

Using an MRAS observer with a dynamic (random) wind speed, the results obtained for Figure 5 
show the evolution of the active power provided by the DFIG with its error as a function of time. Note that 
the curve of the active power P* (in black) followed by the other of the reference (in blue) Pref with a 
transient regime of 0.25 seconds. The maximum peak at 72 KW for a wind speed of 12.2 m/s. It is seen that 
the curve of the active power P* and its error have small fluctuations because of the properties of the 
observer MRAS due to the sensitivity to the dynamic variations of the turbine. 

Under the same conditions described above, Figure 6 shows the evolution of the reactive power 
supplied to the minimum values by the DFIG with its error as a function of time. It is noted that the curve of 
the reactive power Qsmes (in blue) followed the other of the reference (in red) Qsref with a transient regime 
of 0.25 seconds. The maximum value generated from the reactive power is 1.3 Kvar for a wind speed of 12.2 
m/s. This minimum value of the reactive power produced is the consequence of the coupling of the 
mathematical model of the DFIG use. The small fluctuations on the Qsmes curve is the cause of the 
properties of the MRAS observer. 

In the same conditions described previously, Figure 7 presents the dynamic variation of the 
mechanical speed with its error as a function of time. Note that the observed speed curve QMRAS (in blue) 
followed the measured speed QMPPT curve (in red) with a transient period of 0.25 seconds. In this case, the 
transient regime is divided into two parts, the first part where the starting time of the system from 0 to 0.25 
seconds as well as the second part is the delayed time of the observed speed QMRAS with respect to the 
speed measured QQMPPT between 0.25 and 1.5 seconds. This delay is the cause of adaptation mechanism (PI 
regulator) of the observer in which it was favored the accuracy of observation with respect to speed. The 
maximum value of the observed speed is 186 rad / s for a wind speed of 12.2 m/s. The small fluctuations on 
the speed curve QMRAS is the cause of the properties of our observer. 

It can be seen that the reference curves for the active power, which reactivates the mechanical speed, 
has been perfectly followed thanks to the IDA-PBC command and the MPPT command for the dynamic 
variations of the wind. The error curves of the active and reactive powers take important values at startup for 
the transient regime and remain in the average of zero for the permanent regime. 
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5. CONCLUSION 

In this contribution, the MRAS observer and an IDA-PBC controller are developed using advanced 
tools of nonlinear control for wind systems to a DFIG dual feed asynchronous generator. By using both 
MPPT and IDA-PBC commands that guarantee the maximization of the tracked power and the robustness of 
tracking the dynamic variations of the wind. Minimization of tracking errors for active and reactive powers 
are important to ensure reliability of power generation under normal conditions. The results obtained for the 
generated active power and the mechanical speed of the turbine in the dynamic conditions guarantee the 
robustness of our control and the validity of the model proposed, therefore, the precision and the optimization 
to track the variations of the aerodynamic powers. 


APPENDIX 
The parameters used in the model of the DFIG and wind turbine are shown in Table 1: 


Table 1. Parameters used in simulation model. 


DFIG 
Parameter Value Unit 
Pn 100 kW 
Rs 0.455 Q 
Rr 0.19 Q 
Ls 0.07 H 
Lr 0.0213 H 
Lm 0.034 H 
JDFIG 0,53 Kg.m? 
f 0,0024 N.m.s/rad 
p 2 / 
Wind Turbine 
R 35.25 m 
S 1.12 m2 
p 1.22 Kg/m3 
G 90 / 
Jturbine 1000 Kg/m2 
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